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Abstract 

Network theory is rapidly changing our understanding of complex systems, but the relevance 
of topological features for the dynamic behavior of metabolic networks, food webs, production 
systems, information networks, or cascade failures of power grids remains to be explored. Based 
on a simple model of supply networks, we offer an interpretation of instabilities and oscillations 
observed in biological, ecological, economic, and engineering systems. We find that most supply 
networks display damped oscillations, even when their units — and linear chains of these units — 
behave in a non-oscillatory way. Moreover, networks of damped oscillators tend to produce growing 
oscillations. This surprising behavior offers, for example, a new interpretation of business cycles and 
of oscillating or pulsating processes. The network structure of material flows itself turns out to be 
a source of instability, and cyclical variations are an inherent feature of decentralized adjustments. 

PACS numbers: 89.75.Hc,89.65.Gh,47.70.-n,89.40.-a 
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Network theory [1, 2] is an important key to explaining complex systems. Topological 
features [3] determine interesting properties of socio-economic networks [4] or the robustness 
of the world wide web [5], of sensory systems [6], and food webs [7]. Their relevance for 
dynamic features, however, is not well understood. Specific network structures appear to be 
responsible for the beating of a leech's heart [8], pumping effects in some slime molds [9], 
strong variations in ecosystems [10], or glycolytic oscillations in yeast cells [11]. Many of those 
systems, including molecular networks [12], can be viewed as particular supply networks. 
Their functioning is decisive for natural and man-made systems, but prone to oscillations: For 
example, industries often suffer from over-reactions in production [13], disaster management 
struggles with a temporary clustering of forces and materials at some places, while they are 
missing at others [14], and most economies are characterized by "booms" and "recessions" 
[15, 16]. Based on a simple model of interactive flows, we will show that the likely reason 
for these instabilities and oscillations is the underlying network structure of supply systems. 
This is illustrated with a model of macroscopic commodity flows. 

Our simple model of supply networks describes the generation or delivery of products, sub- 
stances, materials, or other resources of kind i at a certain rate Qi(t) > as a function of time 
t. A share X^it) > of that output is required as input for generating or delivering resource 
j, and a share Yi(t) > is absorbed, for example by consumption, degeneration, or loss. If 
Y^j Xij(t)+Yi(t) does not add up to Qi(t), this is assumed to result in a corresponding change 
[Ni(t + At) —Ni(t)] of the stock, inventory, or concentration Ni(t) of resource i per unit time 
At. This assumption implies a conservation equation for resources (which can be general- 
ized). In the limit At — > it reads 



considering that the flow X^t) of resources from supplier % represents a share a^Qjit) > 
of the entire inputs used by supplier j. In Leontief's input-output model of macroeconomics 
[17], the m suppliers correspond to different economic sectors producing one kind of com- 
modity i each. Qi(t) describes the output flow and Ni(t) the inventory of commodity i. The 
stoichiometric coefficients < 1 define an input matrix A = (a^) and reflect the techno- 
logically determined supply network between economic sectors i and j. For simplicity, we 
will treat the input coefficients as constants. 
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The change of the delivery rate Qi(t) in time is generally a function W-i of the stock levels 
or concentrations Nj, their temporal change dNj/dt, and the delivery rates Qj themselves: 
dQi/dt = Wi({Nj}, {dNj/dt}, {Qj}) with j e {1, . . . , m}. In our model of macroeconomic 
output flows, the delivery rates Qi(t) are adjusted in response to two criteria: First, if the 
current inventory iVj(t) exceeds some "optimal" level JV°, the delivery rate is reduced and 
vice versa. A certain inventory level is desireable to cope with variations in the demand. 
Second, if inventories are growing (dNi/dt > 0), i.e. if the current supply Qi(t) exceeds the 
current demand Yi(t) + ^ - a^Qjit), this is an independent reason to reduce the delivery 
rate Qi(t). The following equation for the relative change in the delivery rate expresses these 
responses and guarantees a non-negativity of Qi(t): 

1 dQj = A / iV° _ \_ ji^m 

Qi(t) dt Vl \Ni(t) J Ni(t) dt - U 
Here, z>j is an adaptation rate describing the sensitivity to deviations of the actual inventory 
Ni(t) from the desired one N®. fii is a dimensionless parameter reflecting the responsive- 
ness to relative deviations (dNi/dt) /N^t) from the stationary equilibrium state. Moreover, 
economic systems have an important additional equilibration mechanism, as undesired in- 
ventory levels and inventory changes can be compensated for by adjusting the price Pi(t) of 
commodity %. If the same criteria are applied, we find 

±_dP 1= ( N? \ fM dNj 

Pi(t) dt Vl \Ni(t) J Ni(t) dt ■ U 
We will assume that cej = Vi/vi is the ratio between the adjustment rate of the output flow 
and the adjustment rate of the price in sector %. For simplicity, the same ratio CKj = fii/fii is 
assumed for the responsiveness. 

Increased prices Pi(t) have a negative impact on the consumption rate Yi(t) and vice 
versa, which is described here by a standard demand function fi with a negative derivative 
fi(Pi) = df^/dPf 

Y i (t) = [Y?+amm)). w 

This formula takes into account random fluctuations over time around a certain average 
consumption rate Y® and assumes that the average value of fi(Pi(t)) is normalized to one. 

The possible dynamic behaviors of supply networks can be studied by analytical investi- 
gation of limiting cases (see Appendix) and by means of a linear stability analysis around 
the equilibrium state (in which we have N^t) = Nf, Y^t) = Y?, Q° - Y.j a ijQ°j = Y i \ 
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and Pi(t) = Pf). Supply systems without a mechanism analogous to price adjustment are 
covered by a — > oo and fi(Pi) = const, (no price sensitivity). 

Denoting the m eigenvalues of the input matrix A by uji with \ujj\ < 1, the 3m eigenvalues 
of the linearized model equations are (m times) and 

A i)± » ^ - A ± a/ (A) 2 - 45i ^ , (5) 

where A = ^[Q + a t D t (l - u t )], B % = ^[Q + a t D t (l - oo t )], Q = P?Y?\fl{P?)\/N? and 
Di = QP/Nf. Formula (5) becomes exact when the matrix A is diagonal or the parameters 
fiiCi, oti/iiDi, u,iCi and a^-Di are sector-independent constants, otherwise the eigenvalues 
must be numerically determined. It turns out that the dynamic behavior mainly depends 
on the parameters ctj, z/j//ij 2 , and the eigenvalues u>i of the input matrix A (see Fig. 1): 
In the case ctj — > of fast price adjustment, the eigenvalues Aj^ are given by 2Aj i ± = 
— /ijCj ± a/ (/ijCj) 2 — 4z/jCj, i.e. the network structure does not matter at all. We expect 
an exponential relaxation to the stationary equilibrium for < Vi/fii 2 < Ci/A, otherwise 
damped oscillations. Therefore, immediate price adjustments or similar mechanisms are an 
efficient way to stabilize economic and other supply systems. However, any delay (a, > 0) 
will cause damped or growing oscillations, if complex eigenvalues u>i = Ke(uJi) + ilm(a;j) 
exist. Note that this is the generic case, as typical supply networks in natural and man-made 
systems are characterised by complex eigenvalues (see top of Fig. 1). Damped oscillations 
can be shown to result if all values ^//ij 2 = aiVi/fii 2 lie below the instability lines 

Viln? w {Ci + ttiAfl-Re^i)]} 

{q + q,A[l-Re(^)]} 2 \ 

[aiDihnfa)] 2 J 1 ) 

given by the condition Re(Aj i ±) < 0. For identical parameters Ui/jii 2 = u/ji 2 and ccj = a, 
the minimum of these lines agrees exactly with the numerically obtained curve in Fig. Id. 
Values above this line cause small oscillations to grow over time, but the resulting aggregate 
behavior displays slow variations of small amplitude. The amplitude is limited by phase 
shifts between the oscillations in different sectors and the non-linearities in Eqs. (1) to (4). 
Business cycles, for example, may therefore emerge without exogenous shocks. This result 
markedly differs from the dominating business cycle theory [16]. It is surprising that increas- 
ing oscillation amplitudes are found if the adaptation rates ^ are large. Nevertheless, many 
common production strategies suggest to keep constant inventories A^ , which potentially 
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destabilizes economic systems. Ideal values of z/j//ij 2 should lie below the instability line (6), 
see Fig. Id. 

In some cases, all eigenvalues Ui of the input matrix A are real. This applies to symmetric 
matrices A and matrices equivalent to Jordan normal forms. Hence, the existence of loops 
in supply networks is no sufficient condition for complex eigenvalues Ui (see also Fig. If). 
It is also no necessary condition. Examples for asymmetric matrices with real eigenvalues 
are linear (sequential) supply chains or regular distribution networks. In these cases, Eq. (5) 
predicts a stable, overdamped behavior if all values Ui/ ji 2 = ctiVilfi 2 lie below the lines 

Villi? ~ [Q + aiDi{\ - iOi)]/A (7) 

defined by minj(Aj 2 — 4Bj) > 0. For identical parameters z/«//ij 2 = v/ji 2 and ccj = a, 
the minimum of these lines corresponds exactly to the numerically determined curve in 
Fig. lh. Above it, one observes damped oscillations around the equilibrium state, but growing 
oscillations are not possible. In supply systems without a price adjustment or comparable 
mechanism, Eq. (7) predicts an overdamped behavior for real eigenvalues u>i and V{j ' fi 2 < 
Di(l — Ui) for all i, while Eq. (6) implies the stability condition z>j/ ' fi 2 < D^l — Re(c<jj)]{l + 
[1 — Re(cjj)] 2 /Im(wj) 2 } for all i, given that some eigenvalues uji are complex. 

While previous studies have focussed on the synchronization of oscillators in different 
network topologies [1, 18], we have found that many supply networks display damped oscil- 
lations, even when their units — and linear chains of these units — behave in an overdamped 
way. Furthermore, networks of damped oscillators tend to produce growing ( and mostly asyn- 
chronous) oscillations. Due to the sensitivity of supply systems to their network structure, 
network theory [1, 2, 8] can make useful contributions: On the basis of Eqs. (6) and (7) 
one can design stable, robust, and adaptive supply networks ("network engineering"). For 
example, it is possible to identify structural and control policies which have a dampening 
effect. However, in systems with competing goals (such as intersecting traffic streams), os- 
cillatory solutions can be favourable. The results presented in this study and the applied 
analytical techniques could be also used and generalized to model the dynamics in metabolic 
networks [11, 12] to enhance the robustness of production processes, or to optimize disaster 
management [14]. 



5 



APPENDIX: MATHEMATICAL SUPPLEMENT 



Linearized model equations: The linear stability analysis is based on the following lin- 
earized equations for the deviations rii(t) = N^t) — iV°, pi(t) = P^t) — Pf , and q^t) = 
Qi(t) — Qt from the equilibrium state: 

^ = ?<-E ^ - - > (8) 
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dpi Pf / dn 



% (10) 



eft a? v 

This system of coupled differential equations describes the response of the inventories, prices, 
and production rates to variations £j(i) in the demand. The corresponding eigenvalues are 
shown in Eq. (5). 

Dynamic behavior in limiting cases: Despite their mathematical similarity, Eqs. (2) and 
(3) have a surprisingly different impact on the macroeconomic dynamics: (i) In the case 
oti — > of fast price adjustment, one can eliminate Eq. (10) by assuming q^t) ~ and 
Qi(t) w Q®, so that drii/dt ps Y^\f'{Pf)\pi{t) — Inserting (9) into the time-derivative 
of this equation finally results in the equations 

d 2 n { drii d^ 

—rrr + l^i^i—, h ViCim ~ — (11) 

dt 2 H dt 1 1 1 dt K J 

of damped harmonic oscillators with eigenfrequencies cu® = \Jv{C^ damping constants 7* = 
HiCi/2, and external driving —d^/dt due to variations in the consumption rate, (ii) In the 
case o.i 3> 1 of slow price adjustment or in supply networks for which a price variable is 
not relevant, one can eliminate Eq. (9) by assuming Pi{t) ~ and Pi(t) ~ Pf, so that 
drii/dt ~ qi — a ijQj — Deriving this with respect to time and inserting Eq. (10) delivers 



dt 2 
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where 5ij = 1 for % = j, otherwise 5{j = 0. If we assume sector- independent con- 
stants ctiHiDi = M and a^iDi = V, the 2m eigenvalues \ it ± are given by 2Aj ) ± = 



— M(l — uoi) ± a/[M(1 — uji)] 2 — 4V(1 — uji). For empirical input matrixes A, one never 
finds an overdamped, exponential relaxation to the economic equilibrium, but network- 
induced oscillations (see Fig. 1). An overdamped behavior is only possible if all eigenval- 
ues LOi are real numbers, (iii) If (Ai) 2 / P>i 3> 1, the eigenvalues become Aj 5 _ ps — A { and 
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Aj ;+ ~ —Bi/A,i = —Vi/ni. This corresponds to a relaxation to the equilibrium state in the 
case of a large responsiveness /ii >> 1. An overdamped behavior is found if all eigenval- 
ues uji are real numbers or if all ctj = 0, otherwise one expects network- induced oscilla- 
tions. Interestingly enough, ^ >> 1 implies (dPi/dt)/(niPi) ~ 0, so that Eq. (3) reduces 
to dJVi/cft « ^[A? - Niit)}/^. Therefore, A^(t) « iV° and dQ;/dt « (i.e. Q< « Q?). 
Inserting this into Eq. (1) yields an implicit equation for the price Pi(t) as a function of the 
fluctuations £i(t) in the consumption rate, as usually assumed in economics. It reads 

\y- + tiiiWiim)) ~Qi-J2 = const - ( 13 ) 

j 

Boundary between damped and growing oscillations: Starting with Eq. (5), stability re- 
quires the real parts Re(Aj) of all eigenvalues \ to be non-positive. Therefore, the sta- 
bility boundary is given by maxjRe(Aj) = 0. Writing Cj + aij.Dj(l — Ui) = 9{ + with 
d = P?Y?\f>{P?)\/N» and defining 

6i = Ci + aiDill -Re(^)], 

Pi = ^faiDilm(uJi) (complex conjugate eigenvalues), 

7i = 4z/i//ii 2 , (14) 

we find 



2\/m = -0i - iPt + y/Ri + Hi (15) 

with Ri = 9 2 — j3 2 — ji9i and U = 29iPi — jiPi. The real part of (15) can be calculated via 
the relation 

Re( v / ^±i/ l ) = \J^^/rF+I? + R) ■ (16) 

The condition Re(2Aj///j) = is fulfilled by 7» = and 7, = 46^(1 + 9 2 //3 2 ), i.e. the stable 
regime is given by 

for all i, corresponding to Eq. (6). 

Boundary between damped oscillations and overdamped behavior: For ctj > 0, the imag- 
inary parts of all eigenvalues Aj vanish if lm(ui) = (i.e. $ = 0) and if Ri > 0. This 
requires 

Av- 3 2 

-i = H<9 i -^ = 9 i = C i + aiDi{\ - uj^ (18) 
Hi 2 9i 

for all i, corresponding to Eq. (7). 
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FIG. 1: Properties of our dynamic model of supply networks for a characteristic input matrix 
specified as average input matrix of macroeconomic commodity flows of several countries (top) and 
for a synthetic input matrix generated by random changes of input matrix entries until the number 
of complex eigenvalues was eventually reduced to zero (bottom). Subfigures (a), (e) illustrate 
the color-coded input matrices A, (b), (f) the corresponding network structures, when only the 
strongest links (commodity flows) are shown, (c), (g) the eigenvalues ooi = Re(u>i) +ilm(u;j) of the 
respective input matrix A, and (d), (h) the phase diagrams indicating the stability behavior of the 
model equations (1) to (4) on a double-logarithmic scale as a function of the model parameters 
cti = a and Ui/fii 2 = vjp? = V/M 2 . The other model parameters were set to V{ = C{ = Di = 
Pf = Nf = = 1. Surprisingly, for empirical input matrices A, one never finds an overdamped, 
exponential relaxation to the stationary equilibrium state, but network-induced oscillations due to 
complex eigenvalues Wj. 
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